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Q\ , Abstract 

The fragmentation function of a scalar leptoquark into possible S-wave bound states 
with a heavy anti-quark is calculated to the leading order in perturbative QCD for the 
high energy processes at large transverse momenta. The one-loop equations for the q 2 - 
evolution of moments of the fragmentation function due to the hard gluon emission by 
the leptoquark are derived. The integral probabilities of fragmentation are evaluated. 

> '. 1 Introduction 

An abundant discussion about a possible appearance of new physics at HERA [|T], |3| has 
provided an important and useful experience for a handling of almost-real properties attributed 
^ ! to interactions beyond the Standard Model. So, if a leptoquark [[| is the correct interpretation 

for the excess in the HERA data at high Q 2 , x, then the former has the total width T, which 
^ ' is much less than the QCD-confinement scale, T^q <C Aqcd H- The latter fact means that 

before a decay, the color-triplet leptoquark LQ can be bound with quarks in the lepto-hadrons: 
(gLQ)-baryons or (qiq 2 LQ) -mesons, which are quite exotic states as well as the double LQ- 
^ ! onia: (LQ 1 LQ 2 ). 

It is an attractive and interesting problem to study the spectroscopy, production and decays 
of such hadrons as a possible window of new physics independently of a success or fall off the 
treatment of the HERA events. 

The description of leptoquarks bound with light quarks is a subject of Leptoquark Effective 
Theory, which can be developed as a straightforward continuation of analogous Heavy Quark 
Effective Theory |J with taking a care on the spin structure of leptoquark. The heavy lepto- 
quarkonia: (bLQ) and (cLQ), can be considered in the framework of Non-Relativistic QCD 
keeping in mind again the spin features. 

In this work we discuss the high energy production of heavy lepto-quarkonium containing 
a scalar leptoquark. 

The mo del- independent pair production of free leptoquarks in hadronic collisions was con- 
sidered in H with account for the next-to-leading order QCD corrections. Attaching the result 
to the Tevatron search for the scalar leptoquarks f9j, the authors have found the constraint 
m LQ > 190 GeV. 

At high transverse momenta, the dominant production mechanism for the heavy lepto- 
quarkonium bound states is the leptoquark fragmentation, which can be calculated in pertur- 
bative QCD [pi] after the isolation of soft-binding factor extracted from the non-relativistic 



X 



potential models [O, [L2| . The corresponding fragmentation function is universal for any high 



energy process for the direct production of lepto-quarkonia. 
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In the leading c^-order, the fragmentation function has a scaling form, which is the initial 
one for the perturbative QCD evolution caused by the emission of hard gluons by the lepto- 
quark before the hadronization. The corresponding splitting function differs from that for the 
heavy quark because of the spin structure of gluon coupling to the leptoquark. 

In this work, the LO-scaling fragmentation function is calculated in Section 2. The limit 
of infinitely heavy leptoquark, vtl^q — * oo, is obtained from the full QCD consideration for the 
fragmentation. The splitting kernel of the DGLAP-evolution is derived in Section 3, where 
the one-loop equations of renormalization group for the moments of fragmentation function 
are obtained and solved. The mentioned equations are universal, since they do not depend on 
whether the leptoquark will bound or free at low virtualities, where the perturbative evolution 
stops. The integrated probabilities of leptoquark fragmentation into the heavy lepto-quarkonia 
are evaluated in Section 4 with making the use of non-relativistic wave-functions for the bound 
states. The results are summarized in Conclusion. 



2 Fragmentation function in leading order 

The contribution of fragmentation into the direct production of heavy lepto-quarkonium has 
the form 

da[l H (p)}= / dz da[LQ(p/z),fi] D LQ ^i H (z,fi), 
Jo " 

where da is the differential cross-section of lepto-quarkonium with the 4-momentum p, da is 

that of the hard production of leptoquark with the scaled momentum p/z, and D is interpreted 

as the fragmentation function depending on the fraction of momentum carried out by the 

bound state. The value of /i determines the factorization scale. In accordance with the 

general DGLAP-evolution, the /^-dependent fragmentation function satisfies the equation 

= J z — P L Q^LQ{z/y,fi) D LQ ^i H {y,n), (1) 

where P is the kernel caused by the emission of hard gluons off the leptoquark leg before 
the production of heavy quark pair. Therefore, the initial form of fragmentation function is 
determined by the diagram shown in Fig.[l], and, hence, the corresponding initial factorization 
scale is equal to fi = Ivtlq. Furthermore, this function can be calculated as an expansion in 
a s (2rriQ). The leading order contribution is evaluated in this Section. 




Figure 1: The diagram of leptoquark fragmentation into the heavy lepto-quarkonium. 

Consider the fragmentation diagram in the system, where the momentum of initial lepto- 
quark has the form q = (q , 0, 0, q 3 ) and the lepto-quarkonium one is p, so that 

q 2 = s , p 2 = M 2 . 
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In the static approximation for the bound state of leptoquark and heavy quark, the quark 
mass is expressed as rng = rM, and the leptoquark mass equals m = (1 — r)M. 
The matrix element has the form 



M = -TWm n T ] ™& + ^ - r )lO/V Tfip ~ M)Ih -Mo, (2) 

3V3M 3 r(l — r)(s — m z ) 2 



where the sum over the gluon polarizations is written down in the axial gauge with n 
(1,0,0,-1) 



k ■ n 

with k = q — (1 — r)p. The spinors of Ih and q correspond to the lepto-quarkonium and 
heavy quark associated to the fragmentation. Aio denotes the matrix element for the hard 
production of leptoquark at high energy, R(0) is the radial wave-function at the origin. 
Define 

p ■ n 

z = 

q ■ n 



The fragmentation function is determined by the expression [13 



16tt 2 7 V z 1-zJ |M)| 2 
in the limit of high energies q ■ n — > oo. Then one can straightforwardly find 

D ^ - m iM^r + r2)(1 + (1 " rfz2) ~ 2(1 " r)2(1 + rH (3) 

which tends to 

6fe) = M!Wfez£(4 + M, (4) 
277T uiq r \y° y^J 

at r — > and y — (1 — (1 — r)z)/(rz). The coefficient at the (y — l) 2 term is the same as 
in the fragmentation of heavy quark with the mass m into the S-wave states of the heavy 
quarkonium at y — > 1, if one excepts the factor related with the wave-function of final state. 
The limit of D(y) is in agreement with the general consideration of 1/m-expansion for the 
fragmentation function [H], where 

D{y) = - a (y) + b(y). 
r 

Eq.(f|) determines the a(y) -function explicitly. 



3 Hard gluon emission 



The one-loop contribution of hard gluon emission can be calculated in the way described in 
previous Section. Then the splitting kernel of the leptoquark is equal to 



P, 



LQ^LQ 



>,/i) 



3tt 



■ 2x - 

1 — x_ 



(5) 
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where the "plus" denotes the ordinary action: Jg 1 dxf + (x) ■ g(x) = Jq 1 dxf(x) ■ [g(x) — g(l)]. 
The scalar leptoquark splitting function can be compared with that of the heavy quark 



(x,fi) 



3tt 



1 + x 2 



1 



X 



which has the same normalization factor at x — > 1. 

Further, multiplying the evolution equation by z n and integrating over z, one can get 
from eq. ([3]) the /z-dependence of moments a( n ) for the fragmentation function to the one-loop 
accuracy of renormalization group, 



da 



(n) 



<91n/i 



3tt 



1 



+ . . . + 



1 



n 



1 



n > 1. 



(6) 



At n = the right hand side of (|6|) equals zero, which means that the integral probability of 
leptoquark fragmentation into the heavy lepto-quarkonium is not changed during the evolu- 
tion, and it is determined by the initial fragmentation function calculated perturbatively in 
previous Section. 



a : (n)/a t (2m c ) 
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Figure 2: The QCD-evolution for the averaged fraction of scalar leptoquark momentum, as it 
is developed in the fragmentation with account for the gluon emission to the scale /i, charac- 
terizing the hard production of leptoquark, from the initial value chosen /i = 2m c . 



The solution of eq.(0) has the form 



a (n)(/ i ) = a (n)(/^o) 



Oi s (jj) 



n+l 



.a s (jUo)- 

where one has used the one-loop expression for the QCD coupling constant 

2tt 



(7) 
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Figure 3: The w- function for the leptoquark fragmentation into the heavy lepto-quarkonium 
versus the fraction of r = mg/m. 

where /?o = H — 2n//3 with rif being the number of quark flavors with m q < fi < tjilq- 

Relation (0) is universal one, since it is independent of whether the leptoquark is free 

or bound at the virtualities less than /x . In this work we include the evolution for the 

fragmentation into the heavy lepto-quarkonium. 

As one can see in Fig.§, the leptoquark can lose about 20 % of its momentum before the 

hadronization. 



4 Integral probabilities of fragmentation 

As has been mentioned above the evolution conserves the integral probability of fragmentation, 
which can be calculated explicitly from eq.(H) 



dz D(z) 
w(r) 



8a* \R(0)f 



w(r), 



27 V rriQ 
[7 + 30r + 20r 2 + 20r 3 - 75r 4 



2r 5 + 30r(l + r + 3r 2 + r 3 ) lnrl 



15(1 



(9) 



The function of w(r) is shown in FigFJ at low r. 

To estimate numerically the yield of (JbLQ) and (cLQ) , one has to evaluate the radial wave 
function at the origin from the non-relativistic potential models. As was found in [15, |16| , 
the Martin and Buchmuller-Tye potentials give close results for the heavy quarkonia with the 
reduced mass about m c ~ 1.5 GeV. The characteristics of charmed and beauty lepto-quarkonia 
are presented in Tab.l, where we have used the Martin potential in the limit of infinitely heavy 
leptoquark, m ^> rriq, so that the reduced mass equals the heavy quark mass, and the level 
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energy is given by the sum of quark mass and the binding energy evaluated numerically from 
the Schrodinger equation. 

Table 1: The radial wave-functions at the origin, level energies and average sizes of heavy 
lepto-quarkonia, as evaluated in Martin potential. 



level 


R(0), GeV 3 / 2 


E, GeV 


(r), fm 


lS(cLQ) 


1.61 


1.023 


0.27 


2S(cLQ) 


1.20 


1.606 


0.58 


lS(bLQ) 


3.43 


4.039 


0.16 


2S(bLQ) 


2.56 


4.594 


0.35 



Then one finds for the probabilities of leptoquark fragmentation into the lS-states with 
b and c-quarks: P w (b) = 1.22 • 10~ 4 and P w (c) = 2.16 • 10~ 3 , respectively, at m b = 4.9 GeV, 
m c = 1.5 GeV, a s (2m b ) = 0.18, a s {2m c ) = 0.26 and m = 245 GeV. 

The perturbative fragmentation function in the leading a s -order is shown in Fig.|4] at r = 
0.02. It is quite a hard distribution, which becomes softer with the evolution (see Fig.[|). 

D(z)/N 




1 1 1 1 1 

0.8 0.85 0.9 0.95 1 

Z 



Figure 4: The fragmentation function of leptoquark into the heavy lepto-quarkonium, the iV- 
factor is determined by iV = |^ M jf 2 ^i r )2 , the initial function - solid line, for the (feLQ)-state 
with r = 0.02, the fragmentation function including the evolution - dashed line, at the scale 
a: |Mn^- = 0.25. 
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5 Conclusion 



In this work the dominant mechanism for the production of possible bound states of a scalar 
leptoquark with a heavy anti-quark is considered for high energy processes at large transverse 
momenta, where the fragmentation contributes as the leading term. The corresponding frag- 
mentation function of scalar leptoquark into the heavy lepto-quarkonium can be calculated in 
perturbative QCD, so that for the S-wave states one finds 

fl « - m w& /X~-rw [{l + r2,(1 + (1 " rfz2) ~ 2(1 " r)2(1 + rH 

where r is the ratio of heavy quark mass to the mass of the bound state. In the infinitely 
heavy leptoquark limit, D(z) has the form, which agrees with what expected from the general 
consideration of 1/ m-expansion for the fragmentation functions. 

The hard gluon corrections caused by the splitting of scalar leptoquark are taken into 
account so that the evolution kernel has the form 



4a a (fM) 

Plq->lq{x,II) = — 

07T 



2;r 



1 — X 



which results in the corresponding one-loop equations for the moments of fragmentation func- 
tion (see eqs.(||), (0)). 

The integral probabilities of scalar leptoquark fragmentation into the charmed and beauty 
lepto-quarkonia are of the order of 10~ 3 and 10~ 4 , correspondingly. 
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